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ABSTRACT 

GS 0834 — 430, a 12.3 s accretion-powered pulsar, has been observed in seven outbursts with the 
BATSE large-area detectors on the Compton Gamma Ray Observatory. The first five outbursts observed 
by BATSE occurred at intervals of about 107 days, while the final two outbursts were separated by 
about 140 days. The photon energy spectrum, measured by Earth occultation in the 20-100 keV band, 
can be fitted by a power law with photon index —3.7 or by an exponential spectrum with tem- 
perature kT % 15 keV, with some variations within outbursts. The source has a low pulse fraction, 
<0.15 in the 20-50 keV band. We have observed significant temporal and energy-dependent variations 
in epoch folded pulse profiles. Because the intrinsic torque effects for this system are at least comparable 
to orbital effects, pulse timing analysis did not produce a unique orbital solution. However, confidence 
regions for the orbital elements yielded the following 1 a limits: orbital period P orb — 105.8 ± 0.4 days 
and eccentricity 0.10 < e < 0.17. GS 0834 — 430 is most likely a Be/X-ray binary. 

Subject headings : binaries: general — pulsars: individual (PSR GS 0834-430) — stars: neutron 
X-rays: stars 


1. INTRODUCTION 

GS 0834 — 430 is a transient X-ray pulsar first detected in 
outburst by Granatj WATCH in 1990 February (Lapshov et 
al. 1992a). The source was initially confused with the nearby 
X-ray burster MX 0836-42 (only 24' away and also in 
outburst in 1990 February). Observations by Ginga (Aoki et 
al. 1992), Granatj ART-P (Sunyaev 1991), and ROSAT 
(Belloni et al. 1993) during later outbursts revealed a 12.3 s 
pulse period. The ROSAT observations (Belloni et al. 1993) 
determined the location of GS 0834 — 430 (a = 8 h 35 m 55?l, 
<5 = — 43°11'22", J2000) within a radius of 9" (90% 
confidence). No optical counterpart has been identified. 

WATCH observed six outbursts from GS 0834 — 430 
between 1990 January and 1992 July with peak fluxes of 
300-400 mcrab (8-20 keV). From these data, they derived 
an outburst period of about 1 10 days (Lapshov et al. 1992b). 
The observation of flux from GS 0834 — 430 in the 35-75 
keV band by Granat /SIGMA (Denis et al. 1993), the pres- 
ence of pulsations, and the regular outburst spacing sug- 
gested that this system could be a Be/X-ray binary. 
Although no outburst was observed by Afir/Kvant or 
BATSE at the expected time in 1992 November, A/ir/Kvant 
observed a roughly constant flux of 20-30 mcrab in the 
3-45 keV range from 1992 November 1—16 (Kaniovskii et 
al. 1994). The next outburst occurred about 30 days after 
the predicted time. 
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A total of 10 outbursts from GS 0834 — 430 have been 
observed since 1990 January. Grnnaf/WATCH reported 
three outbursts prior to the launch of the Compton Gamma 
Ray Observatory (CGRO) in 1991 April. Since then, the 
CGRO Burst and Transient Source Experiment (BATSE) 
has observed seven outbursts through continuous monitor- 
ing using Earth occultation and epoch-folding analyses. 
The approximate outburst times determined by BATSE 
detection of pulsations are listed in Table 1. Also listed are 
times spanning observations by other instruments. For pur- 
poses of this paper we will number the outbursts 111, with 
No. 2 denoting the outburst that was not observed. BATSE 
has detected no significant outburst of GS 0834 — 430 since 
1993 June. 

A binary orbital period of 105.8 + 0.4 days has been esti- 
mated from the BATSE timing data. The first five outbursts 
observed by BATSE (Nos. 5-9) occurred at a regular 
spacing of approximately 107 days, in agreement with the 
outburst period derived by the Granat team; however, the 
final two outbursts were spaced by about 140 days and were 
similar in intensity to the earlier outbursts. This behavior is 
unusual since Be/X-ray binaries typically outburst near 
periastron, resulting in regular orbitally modulated out- 
bursts, or they have very large outbursts that are not related 
to the orbital phase (Waters et al. 1989). Interestingly, GS 
0834 — 430 lies below the trend for Be/X-ray binaries on a 
spin versus orbital period diagram (Corbet 1986; Waters & 
van Kerkwijk 1989), in an otherwise “empty” region. 

In this paper, we present the BATSE observations of GS 
0834 — 430. Our observations include a light curve from 
1991 April-1993 June, spectral analyses, and temporal and 
energy-dependent pulse profile variations. We perform a 
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TABLE 1 


Approximate Outburst Times for GS 0834—430 Outbursts Observed by BATSE* 


Start 

Number (MJD) 

End 

(MJD) 

Start 

(Calendar) 

End 

(Calendar) 

Instrument 

Reference 

1 

47908 

47965 

1990 Jan 17 

1990 Mar 15 

Granat 

1 

3 

48164 

48182 

1990 Sep 30 

1990 Oct 18 

Granat 

1 

3 

48208 

48208 

1990 Nov 13 

1990 Nov 13 

ROSAT 

2 

3 b 

48222 

48240 

1990 Nov 27 

1990 Dec 15 

Ginga 

3 

4 

48246 

48293 

1990 Dec 21 

1991 Feb 6 

Granat 

1 

4 

48283 

48285 

1991 Jan 27 

1991 Jan 29 

Ginga 

3 

5 

48374 

48446 

1991 Apr 28 

1991 Jul 9 

BATSE 


5 

48381 

48381 

1991 May 5 

1991 May 5 

ROSAT 

2 

5 

48397 

48404 

1991 May 21 

1991 May 28 

Ginga 

3 

6 

48470 

48542 

1991 Aug 2 

1991 Oct 13 

BATSE 


6 

48592 

48535 

1991 Aug 24 

1991 Oct 6 

Granat 

1 

7 

48590 

48644 

1991 Nov 30 

1992 Jan 23 

BATSE 


7 

48608 

48609 

1991 Dec 18 

1991 Dec 19 

SIGMA 

4 

8 

48698 

48761 

1992 Mar 17 

1992 May 19 

BATSE 


9 

48812 

48842 

1992 Jul 9 

1992 Aug 8 

BATSE 


9 

48818 

48828 

1992 Jul 15 

1992 Jul 25 

Granat 

1 

b 

48927 

48950 

1992 Nov 1 

1992 Nov 24 

Mir / Kvant 

5 

10 

48950 

48998 

1992 Nov 24 

1993 Jan 11 

BATSE 


11 

49101 

49139 

1993 May 24 

1993 Jun 1 

BATSE 



* Times spanning detections by other instruments are also included. 
b Intraoutburst data is included in this interval. 

References. — (1) Lapshov et al. 1992a, 1992b; (2) Belloni et al. 1993; (3) Aoki et al. 1992; (4) 
Denis et al. 1993; (5) Kaniovskii et al. 1994. 


pulse timing analysis and place constraints on the orbital 
parameters. We then discuss the implications of our results. 

2. OBSERVATIONS AND ANALYSES 

BATSE consists of eight identical uncollimated detectors 
positioned on the corners of the CGRO spacecraft such that 
the normal vectors of the detectors are perpendicular to the 
faces of a regular octahedron. The data presented here were 
taken with the BATSE large-area detectors (LADs), which 
are Na I (Tl) scintillation crystals with a geometric area of 
2025 cm 2 and a thickness of 1.27 cm. The LADs are sensi- 
tive to photons from 20 to 1800 keV. The lack of colli- 
mation allows maximum sky coverage for transient and 
burstlike events. Two BATSE data types were used in this 
analysis, the CONT (2.048 s, 16 channel) data and the 
DISCLA (1.024 s, 4 channel) data. The DISCLA data were 
used for pulse arrival time analysis. The CONT data were 
used for spectral fitting, flux and pulse fraction determi- 
nation, and pulse profile analysis. A more complete descrip- 
tion of the instrument and data types can be found in 
Fishman et al. (1989). 

2.1. Source Flux Measurements Using Earth Occultation 

The intensity of a known source is measured by calcu- 
lating the difference in total count rate in source-facing 
detectors just before and after source occultation by Earth. 
Using this technique, BATSE has monitored the flux at the 
location of GS 0834 — 430 continuously since 1991 April. To 
generate the light curve, source count rates were measured 
twice per CGRO orbit using 16 channel 2.048 s time 
resolution CONT data (Harmon et al. 1992). Two minutes 
of data immediately before and after an occultation step 
were fitted with a quadratic background plus a steplike 
function defined by the atmospheric attenuation along the 
line-of-sight to the source. Additional terms in the fit were 
included for occultation steps of known bright sources 
within the fit window. If these bright source steps occurred 
within 10 s of the step being measured, that step was not 


used. Gamma-ray bursts, solar flares, South Atlantic 
Anomaly passages, and electron precipitation events also 
reduced the number of usable steps. Typically, about 10-15 
clean steps per day were available for monitoring GS 
0834 — 430. When the angle between the CGRO orbital 
plane and the source (/? angle) approached 70° (the Earth 
half-angle), the occultations broadened and then ceased 
(Harmon et al. 1992). Every 52 days (CGRO orbital preces- 
sion period), the 28?5 CGRO orbital inclination angle and 
the —43° declination of GS 0834 — 430 combined to a 
maximum greater than 70°, so that the source was not 
occulted for an interval of 3 days. 

Some interfering sources introduced additional com- 
plications to the measurement of GS 0834 — 430. Figure 1 
shows the projections of the limb of Earth at a typical time 
of GS 0834 — 430 rise and set. Detectable sources that are 



Fig. 1. — Nearby sources to GS 0834-430. These sources are variable 
and have occultation times that often occur within the step fitting window. 
The solid curve represents the rising occultation limb + 10 s. The dotted 
curve represents the setting occultation limb ±10 s. The limbs change 
according to spacecraft pointings. At some time during observations of GS 
0834-430, an occultation limb for each of the sources shown was within 
2° of the GS 0834 - 430 limbs. 
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near GS 0834 — 430 are also shown. At intervals equal to the 
spacecraft precession period, the angles between these 
sources and the GS 0834 — 430 occulting limb became small. 
An example of the motion of occultation limbs with time 
can be found in McNamara, Harmon, & Harrison (1995). 
Individual steps where the limb separations were less than 
2 were removed from the flux calculation. This provided a 
more stringent requirement on the data than the 10 s limit 
mentioned earlier. From 1992 August 3-1993 May 22 
(MJD 48840-49100), steps where the limb separation 
between GS 0834 — 430 and the very bright transient, GRO 
J0422 + 32 (Callanan et al. 1995) was less than 5° were 
removed. 

Vela X-l presented an additional problem because it has 
very bright episodes, is only 6 away from GS 0834 — 430, 
and has a 283 s pulse period which is of the same order as 
the step fitting time interval. The long pulse period intro- 
duced steep slopes (which are correlated with occultation 
step size) and relatively short timescale features into the 
background. In an attempt to remove Vela X-l more com- 
pletely, the fit window was extended to 285 s on either side 
of the GS 0834 — 430 step and a Vela X-l template term was 
included, which was obtained from fitting phases in the 
20-50 keV band of BATSE DISCLA data for Vela X-l. 
This method was unsuccessful because it increased the sys- 
tematic errors and residuals from pulse shape variations 
that remained in the data. Most of the effects from Vela X-l 
were removed by averaging standard data over 10 day 
intervals. In addition, outliers that appeared to be associ- 
ated with bright episodes of Vela X-l were removed. 

Occultation measurements from eight energy channels 
(20-160 keV) were averaged over 3 days during GS 
0834 — 430 outbursts and were fitted by a photon power law 
with varying photon index and amplitude. The photon 
index varied from — 3.0 to —4.4 around an average value of 
a % —3.7 and showed no significant trends with intensity. 
The large range of values for a was most likely due to the 
fact that the spectrum was usually constrained by only three 
or four energy channels, so variations in individual channels 
could have affected the spectral fits. The same data were 
also fitted with an exponential model, 


dN 

~dE 



( 1 ) 


with varying temperature kT and amplitude A. The tem- 
perature varied from 12 to 18 keV around an average value 
kT « 15 keV and also showed no significant trends. Again, 
the large range was most likely due to a poorly constrained 
spectrum rather than variations intrinsic to the source. 
Both models were acceptable fits to the data and produced 
equivalent flux estimates. In Figure 2 the Ginga spectral 
model from 1991 May 28 (Aoki et al. 1992) has been over- 
laid with the near simultaneous BATSE data and model. 
The BATSE spectral model is a power law with a photon 
index, ot = — 3.8 + 0.4. The BATSE and Ginga spectral 
models show good agreement. To calculate a bolometric 
correction to the BATSE data, the Ginga model was inte- 
grated from 2-30 keV and the BATSE model was integrated 
from 30-100 keV. The total (2-100 keV) energy flux was 3.5 
times the BATSE energy flux. 

In Figure 3, fluxes were generated by fitting 10 day 
average source count rates with a power law with a fixed 
photon index, a = —3.7. The model was then integrated 



ENERGY (keV) 


Fig. 2. — Combined BATSE and Ginga spectrum for 1991 May 28. The 
dotted line is the Ginga model spectrum (Aoki et al. 1992), and the solid 
line is the BATSE model spectrum. BATSE data points are given by hori- 
zontal bars that are the width of the energy bins. A bolometric correction, 
F total — 3.5F BAXSE , to the BATSE 30-100 keV data was calculated from this 
spectrum. 


from 20-100 keV. The flux values peaked at about 300 
mcrab (where 1 mcrab = 3.1 x 10“ 4 photons cm“ 2 s 1 in 
the 20-100 keV band) and fell to 30-50 mcrab between 
outbursts. Upper limits of 10 mcrab were observed for the 
GS 0834 — 430 location since the outbursts ceased. Hence, 
much of the observed intraoutburst flux was apparently due 
to GS 0834 — 430 plus a small fraction from interfering 
sources. An overall decline in the peak flux can be seen for 
the first five outbursts (Nos. 5-9). However, the next out- 
burst (No. 10), which occurred just after the spacing 
between outbursts increased significantly, was comparable 
in brightness to the first outburst. In addition, a decline in 
the width (as measured by a Gaussian fit to each outburst) 
of the first five outbursts (Nos. 5-9) was observed. 

Systematic errors in the occultation method were esti- 
mated from about 700 days of flux measurements, 1993 
July-1995 May, when GS 0834 — 430 was not detected by 


J u 19 1 Jan92 Jul92 Jan93 Jul93 



Fig. 3. — GS 0834 — 430 flux history in Earth occultation data (20-100 
keV), generated by fitting a power law with a fixed spectral photon index 
a = -3.7 to 10 day average count rates. The error bars shown in the figure 
have been multiplied by 1.65 to include systematic error. The increased 
scatter in the data from days 8840-9100 is caused by flickering of the 
bright transient GRO JQ422 + 32 during outburst (Callanan et al. 1995). 
BATSE has observed no outburst from GS 0834-430 since 1993 June. 
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BATSE. The significances of these fluxes were binned into a 
histogram and the variance, o' 2 , of this distribution was 
calculated. For a Gaussian distribution with no systematic 
error, a = 1. For the 700 days used, a a of 1.65 was 
obtained. The error bars in Figure 3 have been multiplied 
by 1.65 to include systematic errors. 

The flux data, with errors including systematic errors, 
were folded at periods ranging from 1 to 200 days. A j 2 
value for a constant at each period was calculated. In 
Figure 4 a, the period search results for the first five out- 
bursts (Nos. 5-9) are shown. A strong peak is located at 
107.1 ± 0.4 days, and a smaller peak can be seen at half that 
value. However, in Figure 4b , the period search results for 
the last three outbursts (Nos. 9-11), no peak is observed at 
107 days; although a broad enhancement exists near 140 
days and a smaller enhancement occurs at half that value. If 
all seven outbursts are folded, the peak in x 2 is located at 
109.9 ± 0.5 days. The 107.1 day outburst profile for the first 
five outbursts is shown versus arbitrary phase in Figure 4c. 
The outburst profile has no large asymmetries, as also 
observed by Granat (Lapshov et al. 1992b). A power spec- 
trum analysis was also performed on this data set. This 
produced peaks at 105.3 + 0.9 days for the first five out- 
bursts and 114.7 + 1.1 days for all seven outbursts 
(Robinson 1996). 


2.2. Pulsed Observations 
2.2.1. Epoch-folding Analysis and Period History 

An epoch-folding analysis of the 20-50 keV DISCLA 
data was performed for the 800 days containing seven out- 
bursts of GS 0834 — 430. In this analysis, data were rejected 
when GS 0834 — 430 was occulted by Earth, or when solar 
flares, gamma-ray bursts, or other events caused rapid 
variations in the counting rates. The count rates were then 
summed over the detectors exposed to GS 0834 — 430, with 
a weighting for each detector given by the cosine of the 



Fig. 4. — A one-day average flux history was folded at periods ranging 
from 1 to 200 days. A x 2 test was performed, comparing the folded data at 
each period to a constant value. Data folded for the first five outbursts, (a) 
1991 April-1992 August, shows a strong peak at 107.1 + 0.4 days. 
However, data folded for the last three outbursts, ( b ) 1992 July-1993 June, 
shows no strong peaks. ( c ) 107.1 day folded outburst profile for 1991 
April-1992 August. 


aspect angle. Data in 0.5 day intervals were fitted to a model 
in which the pulse source flux was represented by a fifth- 
order harmonic expansion in the pulsar phase. Higher order 
terms were neglected because they could not be measured 
with the 1.024 s resolution DISCLA data. To represent the 
background, the intervals were divided into 300 s segments, 
and the background in each segment was modeled as a 
quadratic in time. The background model value and slope 
were constrained to be continuous at segment boundaries. 
In the calculation of pulse phases from folding the DISCLA 
data, pulse arrival times were corrected to the solar system 
barycenter using the planetary ephemeris Jet Propulsion 
Laboratory DE-200 (Standish et al. 1992), a preliminary 
binary ephemeris (Wilson et al. 1994), and a pulsar rest 
frame phase model. A pulse template was obtained from 
folding the DISCLA data during a strong outburst, subtrac- 
ting the mean, and normalizing to obtain a mean square 
value of 1. The 527 parameterized pulse profiles obtained 
from the 0.5 day interval fits were then correlated with the 
pulse template to obtain a GS 0834 — 430 pulse intensity 
and the phase offset from the preliminary ephemeris for 
each interval. Variations of the pulse shape with time and 
their effects on model fits are discussed in § 2.2.3. 

Figure 5 shows the barycentered pulse period history 
measured by BATSE and the measured periods obtained by 
Granat (Sunyaev 1991 ; Grebenev & Sunyaev 1991 ; Sunyaev 
et al. 1992), Ginga (Aoki et al. 1992), and ROSAT (Belloni et 
al. 1993). Three-day intervals of phase offsets generated 
without assuming a binary ephemeris were fitted by first- 
order polynomials to obtain the periods in the plot. An 
orbital signature is strongly suggested by the presence of 
spin-up and spin-down trends within outbursts. 

2.2.2. Model Fitting and Orbital Constraints 

The phase data were fitted with a binary orbit model plus 
a global polynomial v model. For a first-order polynomial 
in v, a reduced chi-squared value of x 2 = 809 was obtained. 
For higher order polynomials, e.g., seventh-order, x 2 ^ 50. 
Clearly, such a model for the spin-up rate was too simple. 
Next, (p was represented by an orbit plus a different poly- 
nomial in each outburst. For emission times r£ m within the 
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Fig. 5. — Pulse period history of GS 0834-430, corrected to the solar 
system barycenter, with no binary orbital corrections included. Period 
measurements from Granat (Sunyaev 1991; Grebenev & Sunyaev 1991; 
Sunyaev et al. 1992), Ginga (Aoki et al. 1992), and ROSAT (Belloni et al. 
1993) are also shown. 
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ith outburst, 

<r odcl (tn = k + - u + J - o*' . 

(2) 


where i denotes the outburst number, is a pulse emission 
time, T* is a fixed reference epoch in outburst i, and <j) 0i and 
v 0j are adjustable phase and spin frequency at epoch and 


N 

= i 

j = 0 


d u+ i> v (tr - r t .y 

dt^ i) j\ 


( 3 ) 


where d u+ l) v/dt u+ X) are adjustable parameters at epoch t f . 
First, an orbit was fixed and the number of terms increased 
in each outburst until x 2 reached a minimum. Since the 
minimum reduced x 2 = 1.24 at TV = 6, we treated this as a 
measure of noise in excess of counting statistics in the 
phases and increased the phase errors by (1.24) 1/2 to com- 
pensate. Next, an orbit and a different quadratic (JV = 2) 
v model were fitted to each outburst, with all parameters 
allowed to vary. The fit x 2 is written as 


x2 = £ ia - rg (4) 

k = l °<t>k 

for the phase measurements (p k which have errors n 0Jt . The 
phase model </> model is evaluated at the emission time f£ m 
which depends upon both the barycentric observation time 
and the adjustable binary ephemeris parameters. The fit 
converged to a reduced = 1.1 for 487 degrees of freedom, 
but a x sin r, e y and a> were poorly determined. The value a is 
the semimajor axis, i is the angle between the orbital 
angular momentum vector and the line of sight, e is the 
eccentricity, and a> is the longitude of periapse. The param- 
eters of this minimum x 2 value are shown in Table 2. The 
errors on the parameters are 68% single parameter con- 
fidence intervals (Lampton, Margon, & Bowyer 1976). 
Figure 6 shows the measured phases, model components, 
and phase residuals for this model. Figure 7 shows the 
intrinsic v model. One must keep in mind that this is the 
minimum x 2 solution, an example from a large family pos- 
sible solutions. The v model is constrained only by the data. 
However, the presence of a peak in v within outbursts, the 
lack of any global trends, and the concentration of spin-up 
within outbursts strongly suggest that the outburst mecha- 
nism is variable accretion, not obscuration. Variable accre- 
tion is the mechanism believed to power Be/X-ray binary 
systems. 

Since a x sin i and e were not well determined, two- 
parameter confidence regions were examined to constrain 


TABLE 2 

Model Fitting Results at *T n a 


Parameter 

Value 

X 2 (d°0 

534.9(487) 

P „, b days 

105.8 ± 0.4 

T , 

JD 2448810.1!J $ 

a x sin i lt-s 

128!*’ 

ea x sin i lt-s 

15.3!Si 

co° 

140! f| 

/(An 

020!°;?° 


“ Errors are calculated from 68% 
single parameter confidence levels. 



Fig. 6. — Top panel contains the phase points from outbursts Nos. 5-11 
measured by BATSE. The center panel is the two model components of the 
phase model described in § 2.2.2. The dotted line is the intrinsic torque 
component, and the dashed fine is the orbital component. The bottom 
panel shows the phase residuals for this model. 


these parameters. Figure 8 shows confidence regions for a x 
sin i and ea x sin i. The levels shown are 68%, 90%, and 
99%, respectively. Figure 8 shows that although a x sin i is 
not well constrained by the data, ea x sin i is. The pulse 



Fig. 7. — Intrinsic v model evaluated at the minimum x 2 described in 
§ 2.2.2. This is an example of v behavior extracted from a large family of 
allowable solutions. The dotted lines are calculated v values between out- 
bursts. 



0 1 00 200 300 400 

o x sin(i) (lt-s) 

Fig. 8. — Constraints on a x sin i and ea x sin i. Contours shown are 68%, 
90%, and 99% confidence levels. 
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arrival time delay due to a binary orbit is given by (Deeter, 
Boynton, & Pravdo 1981), 

z — a x sin i sin (M -f a>) -f \ea x sin i sin (2 M + co) 

— \e cos co + 0{e 2 ) , (5) 

where M is the mean anomaly. This is used to calculate the 
pulse emission time = t^ SB — z, where f^ SB is the pulse 
arrival time at the solar system barycenter. We note that the 
second term in equation (5) is much better constrained by 
the phase data than the first. 

Next, the values of x p , the epoch of periastron passage, 
and co , the periapse angle, that resulted from the fits within 
the confidence intervals were compared. These two param- 
eters were highly correlated as expected, so a new epoch 
called t 2 was defined. From the second term above, 

2M + « = (t - r 2 ) , (6) 

* orb 

where 


t 2 = 



( 7 ) 


and P orb is the orbital period. Note that t 2 repeats with a 
cycle of P orb or co 4- An. Two-parameter confidence regions 
were examined for t 2 and co . The results of this analysis are 
shown in Figure 9. Again the levels shown are 68%, 90%, 
and 99% confidence, respectively. The plot spans more than 
360° because t 2 has a 4 n cycle. 

Next, we attempted to extract more information about 
the system by replacing a x sin i with the mass function, 


/(Af) = 


An 2 (a x sin i) 

GPL 


03 


( 8 ) 


We then examined two-parameter confidence regions in e 
and /(M). Figure 10 shows that to allow a large mass func- 
tion which would be expected for a Be companion, 
0.04 < e < 0.2, unless the inclination is quite small. If the 
eccentricity is more typical of a Be/X-ray binary 
(0.2 < e < 0.4), then / (Af) < 0.5. 

We have reached several conclusions. The second term in 
equation (5) is well constrained by the data, while the first 
term is not. A circular orbit appears to be ruled out (see Fig. 



0.00 0.10 0.20 0.30 0.40 

Eccentricity 


Fig. 10. — 1 Constraints on /(A/ c ) and e. Contours shown are 68%, 90%, 
and 99% confidence levels. 


8 ). Three of the five orbital parameters (P orb , t 2 , and ea x sin 
i) are well known, but a x sin i and co are poorly known. 

2.2.3. Pulse Profile Variations 

As shown in Figure 11 the GS 0834 — 430 pulse profile 
consists of multiple components which vary considerably 
from one energy band to another. The pulse profiles shown 
are limited to the first three harmonics to minimize noise 
contributions. The error bars are spaced such that individ- 
ual points are uncorrelated. The leading component 
(^ phase 0.3-0.75) appears to be dominant or at least equiv- 
alent in the 20-30 keY band, while the trailing component 
(^ phase 0.75-0.3) is dominant at higher energies. Pulse 
profiles observed by BATSE are similar to those observed 
by Ginga in the 1.2-27.5 keV range (Aoki et al. 1992). Both 
BATSE and Ginga observed significant pulse shape varia- 
tions with energy. The pulse profiles from different time 
intervals are not phase aligned. Hence any shift in phase 
between the different time intervals shown in Figure 1 1 is 
arbitrary; however, the apparent shift in phase of profile 
features between different energy bands within a single time 
interval is intrinsic to the source. 

CONT data were folded to obtain pulse profiles. The 
phase model used consisted of a quadratic and an orbital 
model obtained from DISCLA data. Due to the large intrin- 
sic v, which was not simply behaved, only short ( <2 weeks) 



Fig. 9. — Constraints on x 2 and (o. Contours shown are 68%, 90%, and 
99% confidence levels. 



RELATIVE PHASE 

Fig. 11. — ( a ) Pulse profile from outburst No. 6, 1991 September 19-23 
(MJD 48518-48522), used as the first template in temporal pulse shape 
variations analysis, (b) Pulse profile from outburst No. 7, 199i December 
28-1992 January 2 (MJD 48618-48622). 
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Fig. 12.— Pulse fraction calculations for GS 0834 — 430, using 
(r pu i s « d - ^ u »«d m J/'occ’ where r is the count rate. 


intervals were used to minimize smearing in individual 
pulse profiles. Any trends remaining in the phase residuals 
were removed by fitting 1 day bins with a quadratic which 
was then used to correct the pulse profiles. 

The algorithm used to compute pulse profiles introduced 
correlations between adjacent phase bins, making a sta- 
tistical comparison of different pulse profiles difficult to cor- 
rectly formulate in the time domain (see Appendix). 
However, such a comparison is easily formulated in the 
frequency domain. As shown in the Appendix, the harmonic 
coefficients of the pulse profile are uncorrelated (although 
their errors are dependent upon harmonic number). These 
harmonic coefficients were calculated as 

a k = Tj £ Rje- 2nijkl \ (9) 

/V j=l 

where i = (— 1) 1/2 , TV = 64 is the number of phase bins, and 
Rj is the count rate in the jth bin of the pulse profile. Only 
the first three harmonics were kept because higher order 
harmonics were not statistically significant. The variance on 
the harmonic coefficients of order k is 


ffRe(«)=y 0°) 

where g\. is the variance of R } and b k is a correction due to 
the rebinning. A detailed derivation of b k is given in the 
Appendix. 

A study of pulse profiles from all seven outbursts 
observed by BATSE has resulted in the discovery of two 
different pulse shapes. For this analysis, harmonic represen- 
tations of 3-6 day summed pulse profiles were compared. 
The time interval lengths were chosen to maximize signal, 
while minimizing phase drift. Only the first three harmonics 
were considered in the analysis because the higher harmo- 
nics were not statistically significant. A template was chosen 
from the set of pulse profiles, and it was compared to all 
available profiles by minimizing 


,2 _ 


* * \aj-t c k p c e 2 ^ k \ 2 

c=l k = 1 l>R< (a k c )] 2 + ft 2 K c((1 «)] 2 


(ID 


where k is the harmonic number, c is the energy channel 
number, is the harmonic coefficient of the template profile 
for channel c, a c k is the harmonic coefficient of the compared 
profile for channel c\ $ is a phase factor common to all 


energy channels, and is an energy channel-dependent 
scaling factor. The time interval from 1991 September 
19-23 (Fig. 11a) was chosen as the initial template and 
compared in each energy channel with 40 intervals span- 
ning the seven outbursts. Intervals with significant signals in 
less than three energy channels were not included in the 
comparisons. Of the 40 intervals studied, three showed sig- 
nificant deviations from the template at the 99% confidence 
level. Next, the most significantly different interval, 1991 
December 28-1992 January 2, shown in Figure 116, was 
chosen as the new template. We found 13 intervals that 
were different from it at the 99% confidence level. No inter- 
vals differed significantly from both templates. The tempo- 
ral pulse shape variations were not clearly correlated with 
intensity or orbital phase. 

In the phase analysis described in § 2.2.1, the same tem- 
plate was used for all outbursts. The template used was 
generated using DISCLA data which corresponded approx- 
imately to combining the first three energy channels in 
Figure 11a. To estimate the effects of pulse shape variations 
on the orbital model, we first compared the pulse period of 
12.3 s to a x sin i = 128 lt-s. The pulse period is approx- 
imately 10% of a x sin /. Therefore, if the pulse shape varia- 
tions produced a 0.25 cycle shift in the measured phase 
(which we think is unlikely), the effect on a x sin i would be 
only about 2% which is well within our contours. In a 
system as wide as this one, the effects of pulse shape varia- 
tions on the orbital solution are fairly small. 

Pulse fraction calculations were performed by calculating 
the mean minus the minimum value for the epoch-folded 
pulse profiles and dividing that by the average occultation 
count rate for the interval. Figure 12 shows pulse fractions 
calculated for five different intervals from seven outbursts. 
Clearly, the pulse fractions are quite small, <0.15 for 20 50 
keV. Also, the pulse fractions show a marginal increase with 
energy. 

3. DISCUSSION 

The most likely scenario describing GS 0834 — 430 is a 
Be/X-ray binary with a wide orbit, a small (but nonzero) 
eccentricity, and a low inclination angle. Figures 8 and 10 
rule out a circular orbit. A low inclination angle would give 
a small mass function and a low pulse fraction (Fig. 12), 
assuming the spin and orbital angular momentum vectors 
are aligned. These are both consistent with observations for 
GS 0834 — 430. For a typical Be star of 15 M e , an inclina- 
tion angle of 20° gives a mass function of 0.5, which is 
within the 68% contour in Figure 10. The first nine out- 
bursts were periodic and similar in intensity, consistent with 
normal (class I; Stella et al. 1986) outbursts observed from 
other Be/X-ray binaries. In addition, the peak flux (Fig. 3) 
shows a systematic decline for outbursts 5 9. The last two 
outbursts (10 and 11) were shifted in orbital phase, with the 
largest shift occurring for the final outburst. Also, these 
outbursts did not show the declining intensity trend of the 
previous five. Thus, the first nine outbursts appear to be 
related and separate from the last two outbursts. A change 
in the environment around the neutron star that is not 
understood apparently occurred between outbursts 9 and 
10. GS 0834 — 430 is unusual because the combination of 
periodic and aperiodic outbursts with similar intensities has 
not been observed in other systems. The periodic transient 
behavior is similar to some Be/X-ray binaries, although no 
companion has been found. Optical CCD images showed 
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three stars lying within the ROSAT 9" error circle with 
R < 23.5 (Belloni et al. 1993). Only the brightest star was 
studied further, and its spectrum is typical of a normal late- 
type star. Additional optical and infrared studies are needed 
to confirm the hypothesis that the companion is a Be star. 

Observations of steady spin-up during outbursts allow us 
to calculate a lower limit to the distance to GS 0834 — 430 
(Chakrabarty et al. 1993). The observed values of the pulsar 
spin period and derivatives during periods of maximum 
spin-up imply that M < 1.1 x 10 17 g s -1 . Combining this 
with the concurrent 30-100 keV flux observations and the 
bolometric correction F total = 3.5F BATSE calculated in § 2.1, 
we infer that the distance to GS 0834 — 430 is most likely 
>4.5 kpc. 

We close by mentioning two other X-ray binary systems 
that exhibit characteristics similar to GS 0834 — 430. The 
first of these systems provides a possible scenario which 
would allow a Be companion in a nearly circular orbit. The 
second system suggests possible transient outburst mecha- 
nisms that could allow relatively periodic outbursts that are 
not dependent on periastron passage. For 2S 1553 — 542, a 
system generally accepted as a Be/X-ray binary with a 
nearly circular orbit, Kelley, Rappaport, & Ayasli (1983) 
suggest that a small but finite eccentricity could be produc- 
ed by a supernova explosion involving a small mass loss 
(<2M 0 ) from a less than 4 M 0 He star. This model places 
a mass limit of less than 20 M 0 on the companion’s present 
mass, allowing a Be companion. GX 301—2, which is a 


high-mass supergiant system in an eccentric orbit, shows 
regular near-periastron outbursts and transient weaker 
apastron outbursts (Pravdo et al. 1995; Koh et al. 1996). 
Pravdo et al. (1995) suggest that the outbursts occur 
because the orbital plane of the neutron star and the equa- 
torial disk of material surrounding the companion are mis- 
aligned. Perhaps GS 0834 — 430 outbursts could be 
described by this model, with a suppression of the second- 
ary outbursts early on followed by an enhancement of these 
secondary outbursts. Haberl (1991) gives a different possible 
model for GX 301—2 which consists of a gas stream caused 
by a region of enhanced mass loss toward the neutron star 
on the surface of the primary. Outbursts occur when the 
neutron star passes through the spiral shaped gas stream. 
Different orbital geometries could change the phases of the 
outbursts. Perhaps an instability in the gas stream could 
cause behavior like that seen in the last two outbursts of GS 
0834 — 430. Future observations of GS 0834 — 430, if it reap- 
pears, could improve the orbital parameters and possibly 
give insight into the behavior of other systems. 
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5-3109. We thank an anonymous referee for helpful com- 
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APPENDIX A 


APPENDIX CORRELATIONS DUE TO REBINNING 


A pulse profile is easily constructed from a series of photon events by assigning each arrival time t k a pulse phase according 
to <j>(t k ) = vt k mod 1, where v is the pulse frequency, and by building up an N - bin histogram of phases Sj where j = 1, . . . , N. If 
the time series is itself binned, one must account for the fact that individual time series bins may overlap more than one pulse 
profile phase bin. One approach is to treat each time series bin as a delta function at the bin center, and add it entirely into the 
single appropriate phase bin (whole bin folding). While this has the advantage of maintaining statistical independence of the 
phase bins, it results in some loss of phase information, especially if l/2Vv < t, where z = t k + l — t k is the time series spacing. 
An alternative approach is to split an overlapping time series bin over each of the overlapping phase bins in proportion to the 
degree of overlap, making the assumption that the accumulated counts were uniformly distributed over the time series bin 
( fractional bin folding). While this will maximally retain phase information, it introduces correlations between phase bins, 
resulting in apparent fluctuations which can falsely mimic a pulsed signal if not properly accounted for. In this Appendix, we 
compute the statistical correlations introduced by rebinning. 

Consider a binned time series with C k counts in the /cth bin at midtime t kf with variances a k and uniform time spacing t. We 
can epoch fold this time series into an N - bin pulse profile with midbin phases $ 7 - = (/ + j)/N by computing each phase bin as 

S J = I fjk Ck , (Al) 

k 

where the rebinning factor/}* for fractional bin folding is given by 

j r*k+ 1/2 

fjk = -\ , (A2) 

T Jtk-z/2 

with 


Fj{(t>) = 


1 for (pj 
1 0 otherwise , 


A<b Acb 

Y < m <<t>j + Y 


(A3) 


where A<f> = \/N. {Note that for whole bin folding, f jk = F,[0(t ft )].} Similarly, the integration time accumulated in each phase 
bin is 



(A4) 
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Fig. 1 3.— Calculated and Monte Carlo simulated values of the rebinning factor b k used to account for correlations in pulse profiles. The Appendix 
describes b k in detail. 


o Calculation 
o Monte-Carlo 


so that the count rates in the pulse profile are just Rj = Sj/Tj . The square pulse function F can be represented by the 
infinite Fourier series, 

FJL4>) = jj t sine , (A5) 

where sine x = sin x/x and i = (— 1) 1/2 . Since the instantaneous pulse phase is = (p(tj) + v(t — t } ), where v is the pulsar 
spin frequency, the rebinning factor can be rewritten as 

f jk = — Y s * nc (rcmvt) sine \ . (A6) 

The covariance of the pulse profile rates is 

Cov (Rj, R k ) = X ’ < A7 ) 

n l j l k 

where f )n and f kn are computed according to equation (A6). If we assume <f>(t k ) is uniformly distributed (which is equivalent to 
requiring that the length of the time series is > 1/v and that x and the pulse period are incommensurate), then the phase- 
averaged covariance is 

Cov ( Rj , R k ) = X sine 2 (mnv t) sine 2 ( • ( A8 ) 

N Tj T k m = _ x \ N / 

We can rewrite this in terms of a correlation coefficient, 

Cov (Rj, R k ) _ sine 2 (mnv x) sine 2 (mn/N)e 27zim[ * (tk) -M 

[Var (Rj) Var (R*)] 1/2 Ym = -<» sine 2 (w7tvt) sine 2 (ntn/N) 

where Var (Rj) — Cov (R j9 Rj) is the variance of an individual pulse profile bin. 

We are now in a position to examine the statistical behavior of the pulse profile in the frequency domain by means of 
harmonic decomposition of the pulse. The harmonic coefficients of the pulse profile are 

= i I Rje- 2 ’**', (A 10) 

iv J=1 

and the covariance between these harmonic coefficients can be calculated as 


Cov (a k , a,) = £ I Cox (Rj, R m )e- 2 "« k *>- 1 *”' . 

iV j=l m=\ 


Substituting for Cov (R jy R m ) using equation (A9), we find 


Cov (a k , a,) = S kl Var ( R k ) 


sine 2 l(k + ;W)7rvt] sine 2 [( k + jN)n/N~\ 
£*= - <x> sine 2 (m7tvt) sine 2 (rin/N) 


from which it follows that the harmonic coefficients are uncorrelated and that 


<7Re(a t ) = = 1 Var (d k ) . 


(All) 


(A 12) 
(A 13) 
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We can write the variance of the real part of harmonic coefficient as 

^ Var(R) , (A 14) 

where the rebinning correction factor b k is given by 

_ N IT— Sine 2 [(fc + j'N)7tvt] sine 2 [(fc +;/V)7t//V] 

* „ sine 2 (nnvT) sine 2 (rtn/N) 

Note that for whole bin folding, b k = 1. For the fractional bin folding of CONT data used to study GS 0834 — 430 in this 
paper, we have N = 64 phase bins and x = 2.048 s, and 1 / v = 12.3 s. This calculation has been verified by Monte Carlo 
simulations using these values. The results of the calculation and simulations are shown in Figure 13. 
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